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Abstract 

We analyze multiparticle production in a thermal framework for 7 cen- 
tral nucleus nucleus collisions, e + + e~ annihilation into hadrons on the 
Z resonance and 4 hadronic reactions (p+p and p+p with partial central- 
ity selection), with center of mass energies ranging from \fs= 2.6 GeV 
(per nucleon pair) to 1.8 TeV. Thermodynamic parameters at chemical 
freeze-out (temperature and baryon and strangeness fugacities) are ob- 
tained from appropriate fits, generally improving in quality for reactions 
subjected to centrality cuts. All systems with nonvanishing fugacities 
are extrapolated along trajectories of equal energy density, density and 
entropy density to zero fugacities. The so obtained temperatures ex- 
trapolated to zero fugacities as a function of initial energy density e i n 
universally show a strong rise followed by a saturating limit of Tn m = 
155 ± 6 db 20 MeV. We interpret this behaviour as mapping out the 
boundary between quark gluon plasma and hadronic phases. The ratio 
of strange antiquarks to light ones as a function of the initial energy 
density e j n shows the same behaviour as the temperature, saturating 
at a value of 0.365 ± 0.033 ± 0.07. No distinctive feature of 'strangeness 
enhancement' is seen for heavy ion collisions relative to hadronic and 
leptonic reactions, when compared at the same initial energy density. 
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1 Introduction 



Hadronic reactions involving copious production of secondary particles have 
been associated with an underlying thermodynamic behaviour since the earliest 
observations in cosmic rays |l]]. The observable energy regime ranging up to 
140 TeV per nucleon pair corresponding to incident cosmic ray nucleons of 
E < 10 10 GeV demands elaborate simulations of the observed extended air 
shower development , in order to extract definite multiplicity distributions of 
the elementary hadronic and nuclear reactions ||. 

Thermodynamic models are widely and successfully used to describe identified 
particle ratios in hadronic and especially heavy ion collisions |J , [[|, [[J, 0, 
0, but were also extended to e + e~ annihilation into hadronic final states || 
, 0. In heavy ion collisions two main parameters - energy per nucleon pair 
and centrality - are varied and their influence on thermodynamic variables - 
temperature and chemical potentials - are studied. 

We derive and discuss thermal properties in search of the phase boundary 
between quark gluon plasma and condensed hadrons. We consider as one 
characteristic signature of this boundary the critical dependence of kaon num- 
ber densities on the initial energy density. The kaon multiplicities observed 
in Pb+Pb collisions at y/s=17 GeV [[UJ, and in other nucleus+nucleus and 
p+p collisions at y/s ~ 5-19 GeV |TT|JT^| , serve as motivation for our thermal 
description extending the previous analysis of one of us (SK) |Tl|Jl^| (see figure 

ID- 

In section 2 we discuss the dependence of the critical temperature on the vac- 
cum pressure in QCD, limited to the case of vanishing chemical potentials, 



extending the work of reference ||13| . In section 3 we seek to assign each ther- 
modynamic state established at chemical freeze-out for any given reaction for 
which chemical potentials for baryon number and strangeness do not vanish, 
an equivalent state at zero fugacities. We do this extrapolating along curves of 
equal entropy density, energy density and number density. The extrapolation 
of the temperature of systems with finite fugacities to zero fugacities, and the 
investigation of their dependence on the initial energy density has been pro- 



posed in [ 12| . We find that a more universal parameter as a measure of positive 



and negative strangeness production is A s = 2 ^/-7\ ■ This parameter is 
also extrapolated to the zero fugacity systems. 

The details of this procedure are worked out in the first four subsections of 
section 3. The results at zero fugacities are contained in subsection 3.5 and 
represented in figures ||| - 1|] in relation with the energy density initially achieved 
in each reaction. Despite large errors, the phase boundary between quark gluon 
plasma and hadronic phases is clearly mapped out. 
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We use in the hadronic phase the approximation of noninteracting hadron 
resonances to describe in this sense global ratios of hadrons produced in the 
following reactions : 



1) central Au+Au collisions at RHIC ^ = 130 A GeV , p| . 



2) central Pb+Pb collisions at SPS at Vs = 17 A GeV p 



central S+A collisions ( A=Au,W,Pb) at SPS y/s = 19 A GeV [T7| 
central S+S collisions at SPS = 19 A GeV p] . 



3) central Si+Au collisions at AGS y/s = 5.4 A GeV | 



central Au+Au collisions at AGS xfs = 4.9 A GeV pf, J2G 



4) Ni+Ni collisions at GSI yfs = 2.8 A GeV MM ■ 

5) e + + e- collisions at LEP y/Ts =91.19 GeV g, § 



6) p+p collisions at t/s = 27 GeV [^3| , 
p+p collisions at y/s = 17 GeV pi , p 



7) p + p collisions at yfs = 900 GeV 



peripheral p + p collisions at \/s = 1.8 TeV p 



central p + p collisions at \fs = 1.8 TeV |26 . 
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2 Outline of basic points 



We consider first the thermodynamic potential $ of a grand canonical en- 
semble of hadron resonances without further interactions among them. The 
thermodynamic variables are 

V, T, x b = V b I T, x s = V s /T (1) 

i.e. volume, temperature, baryon and strangeness fugacity. The fugacity of 
the third component of isospin is set to zero, neglecting isospin asymmetries. 

$ = gV;g = g ( (3 , x b , X s) = E a 9 a + 9 o 

(2) 

(3 — 1 / T ; g = (3 p ; p : pressure 

In eq. @ the sum extends over hadron resonances denoted by a, where we 
only include the pseudoscalar and vector u,d,s meson nonets as well as the spin 
1/2 baryon octet and spin 3/2 decuplet and their antiparticles as well as the 
/ (400 — 1200) or a, interpreted as scalar glueball |27J for simplicity. 



go — (3 P o takes into account the nonzero vacuum pressure characterizing 
the hadronic phase in QCD, which we restrict to the three light flavors. 



9 

Po (T = 0) = B 2 - \K{m q 

32 7T 2 



0.00302 GeV 4 
0.00658 GeV 4 



B 2 = ( | \ F^ v F a ^ | ) 



0.125 GeV 4 [EB1 



0.250 GeV 4 |29] 



A = E g ™</(0|-gg|0)~ f l(\ml + m 2 K ) = 0.00217 GeV 4 

(3) 

Converting energy density units to GeV / fm 3 the range of values for the 
(positive) vacuum pressure becomes 



Po (T = 



0.377 GeV / fm 3 

0.823 GeV / fm 3 (4) 
PoCO ~ PoCT = 0) (1 - (T/T cr ) 4 ) 
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The temperature variation of the vacuum pressure is an approximation with 

T C r -^cr ( X 6 i X s )• 

The quantities g a in eq. (fj) for noninteracting resonances are then given by 

/oo 
dE E p I ( 2vr 2 ) / 
m a 

I = T log [ 1 T exp ( - f3 E + x 6 B a + x . S a ) ] ( 5 ) 
w a = (2/ Q + l)(2Sp Q + l);p = 2 _ m 2 

In eq. (|j) ( / , , 5 , S ) a denote isospin, spin, baryon number and 
strangeness of the hadron resonance a. The -,+ signs apply to bosons and 
fermions respectively. 

If we take into account the variation of the vacuum pressure with temperature 
then the masses m a — > m a ( T ) of quasiexcitations become temperature 
dependent quantities. 

In all subsequent calculations we neglect these temperature dependent effects, 
which however do set in dramatically, when the temperature deviates from the 
crirtical one by less than 10 %. 

p o ( T ) -> p o = p o (T = 0) ; m a (T) -> m a =m a ( T = ) 

(6) 

We state here, that a thermal description cannot account for any azimuthal 
( (p - ) dependence of inclusive cross sections. Thus azimuthal anisotropies 
provide a measure for non-equilibration. 

The potential $ in eq. (0) and the entropy denoted by 5 give rise to the 
differentials 

d $ = - S d (3 + N v dx v + 9 dV 
dS = + (3 dE - x v d N v + g dV -> 

(7) 

< s = ^ s y = $ + /3£-x^^^; is = b , s 

Q s = (3(e + p)-XuQu 

In eq. (0) s denotes the entropy density, e the energy density and g b, s the 
net baryon and strangeness (number) density respectively. 
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Qv = J2 a V a Qa;Va = (B a ,S a ) (8) 

As is apparent from eq. ([7]) the entropy density obtains no contribution from 
vacuum energy density and pressure, as long as the vacuum pressure is taken 
as independent of temparature as indicated in eq. (|6]). 
Energy density and total number density then take the form 



= Q e ~ P ; e = E a Q ea 5 Q n = E a Q 

Q ea I f 00 | E 

w a dE E p / {2n 2 ) I \ n 

0<>» I J ma i 1 I ( 9 ) 



n = 

exp ( p E - x b B a - x s S a ) =F 1 

The approximation of hadronic interactions by free hadron resonances is of 
course at best a consistent approximate thermodynamic description. We will 
use it here to describe in this sense global ratios of hadrons produced in the 
reactions listed in the introduction. 

We will present a detailed analysis of these collisions, as seen at the time of 
chemical freeze-out in the next section. The key parameter characterizing not 
this freeze-out instance but the time of first onset of thermal conditions is the 
initial energy density £j n . It determines the thermal parameters of the hadronic 
system as seen at chemical freeze-out and distinguishes those systems which 
have hadronized after transcurring the quark gluon plasma phase from those 
which remained throughout in the hadronic phase. 
The phase boundary for vanishing fugacities 

In the remainder of this section we analyze the equilibrium condition for both 
phases to coexist in the case of vanishing fugacities or equivalently chemical 
potentials. This latter choice is just for clarity of argument and also because 
for limiting center of mass energies and finite baryon number, strangeness and 
charge of the initial system this is the limiting case. 

This condition corresponds to equating the Gibbs densities of the two phases 
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H 



T 



9 QGP (T); Xu 







T 



T , 



Z a g a + Ppo ; g QGP = Z- a g-« 

+ : hadron resonances 
— : gluons, quarks and antiquarks (u,d,s) 



(10) 



In the same approximative vein as applied to the hadronic phase (H) we ap- 
proximate the various quasiexcitations in the quark gluon phase by free gluons, 
quarks and antiquarks analogous to the form of g a in eq. (H) 



(j -,, '■" -„ / dE E p I ( 27r 2 ) I 

% —ot 

I = T log [ 1 T exp ( - E ) ] ; p = ^E^ 

gl : w g i = 16 m g i = - 
q : w q = 6 m g + 
g:w 5 =6 m g + 



m 



-a 



(11) 



5.25 , 8.75 , 175 ) MeV 



In eq. ( |TTD the color/spin weights, masses and boson/fermion signs are spec- 
ified. The quark masses represent a renormalization group invariant 'best' 
choice. It turns out that near the critical temperature T cr = O ( 200 MeV ) 
the (anti) quark contributions to the Gibbs density do not deviate substantially 
from the limit of vanishing quark masses including m s . 

The crossings of g H and g QGP as displayed in figure [I] seem to suggest a first 
order phase transition with respect to energy density. 

This is however due to the approximation of fixed masses for hadron excitations 
and of free quark and gluon flavors in the region of T cr = 194 ± 18 MeV. 
We note that the above estimate of the critical tempeature at zero fugacities 



confirms within the theoretical error the previous estimate of one of us [13 



Within the approximations and errors we obtain the critical temperature and 
critical hadronic energy density on the side of the hadronic phase for vanishing 
fugacities 
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C5 
O 
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10 



0.08 



'B2=0.125_GeVM' 
'q_qbar_g' 
■B2=0.1875_GeVM' 
'B2=0.25_GeVM' 




0.1 



0.12 



0.14 0.16 0.1E 
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0.2 0.22 0.24 0.26 



Figure 1: Gibbs potentials as a function of the temperature for three values 
of the gluon condensate in the ground state, and for free quark flavors u, d, s 
and gluons. 



T cr = 194 ± 18 MeV 



I 212 MeV \ 

194 MeV 
\ 176 MeV ) 



1.797 

1.018 GeV/fm 3 
0.5406 

(12) 



Obviously the critical boundary curve could be extended to arbitrary chemi- 
cal potentials, but again the 'melting' of masses of hadronic quasiexcitations 
modifies the Gibbs density on the side of the hadronic phase, whereas the ap- 
proximation of free quark and gluon modes does not warrant the extrapolation 
to arbitrary chemical potentials, especially for small temperature. 
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3 Reduction of chemical freeze-out parame- 
ters to zero fugacities in multiparticle pro- 
duction 

In this section we apply the thermal model introduced in the last section 
(equation 5) to extract the intensive thermodynamic parameters defined in eq. 
1. We then extrapolate these parameters to zero fugacities along states with 
equal entropy density, energy density or density The goal of this analysis is 
to compare the temperature at zero fugacities with the initial energy density 
achieved in the collision in order to reveal a boundary reflecting the QCD 
phase transition. 

One serious well known problem when comparing models to experimental data 
is due to decays of resonances. This effect is called 'feeding'. We compare cal- 
culated particle ratios to experimental data, taking all (strong, electromagnetic 
and weak) decays with a branching ratio not below 1% into account. We try to 
account for experimental acceptance for K + and K~ , assuming a 50% feeding 
to pions, due to their long decay length. The other weak decays (e.g. A, K®) 
have a much shorter decay length and are assumed to fully feed into secondary 
particles. The K? decay is not considered. 



Contrary to references we do not take all hadronic resonances below 2 
GeV mass into account. The theoretical error inherent to the free resonance 
approximation allows in our opinion our reduced set. This is supported by the 
finding that we derive similar thermodyamic parameters with similar accuracy 
as e.g. the ones found in reference [TjJ. Though the high mass resonances 



are definitely produced at sufficiently high energy, it is not clear that the 
thermal description becomes more accurate by including them because the 
quasiexcitations may not correspond to them, especially if these are quark and 
gluon modes. All similar thermal descriptions suffer from replacing interactions 
by noninteracting resonances. Furthermore, the isotropic angular distributions 
characterizing the models are far from reality 0J31|] (e.g. flow phenomena, 
diffractive phenomena). 

We enforce strangeness conservation, to obtain the best fit, exept in the case of 
p+p collisions at y/s—17 GeV. We find the error on the fitted parameters, at 
± 1 unit of the ratio of y 2 / DOF away from the best fitted value. Our quoted 
error reflects the inherent theoretical error, despite the difficulty to quantify 
it. Note that other analyses finding a smaller error seem to use as error 
the fitted parameters at ± 1 unit of the x 2 itself. 
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In the following sections we note if the data we use are taken with a minimum 
bias trigger (that means no trigger bias is imposed) or a central respectively a 
peripheral trigger. Central is here understood as the lowest impact parameter 
region. It is usually selected in the experiment by considering only collisions 
with the largest particle multiplicity (e.g. CERN-WA97, Fermilab-E735) or 
the largest transverse energy (e.g. CERN-NA52), or the smallest forward go- 
ing energy (which reflects the spectator nucleon number) (e.g. CERN-NA49). 
Respectively, peripheral collisions are a selection of collisions with the largest 
impact parameters using the same means as described above (charged multi- 
plicity, energy). 

In the following analysis we also give the quantity: A s = ^1 which we take 
as a measure of the strangeness suppression factor defined and used in the 
literature as X s = ^^jj^ 0- We use antiquarks to consider in a simplified 
way only the newly produced valence antiquarks. We give also the equivalent 
\ s at zero fugacities. From quark counting rules one expects A s to be in the 
range 0.3-0.5 In principle one should include in any definition of A s the 



newly produced sea quarks however the relevant proportion of the latter inside 
hadrons, is still subject of experimental and theoretical investigations. 

Generally we find that thermal fits to minimum bias or peripheral nuclear and 
hadronic reactions are not satisfactory, in contrast to central collisions. We 
conjecture that this is due to the presence of at least two thermal sources which 
we attribute to diffractive versus pomeron induced subprocesses. Diffractive 
processes at finite \fs feed back into the midrapidity region, whereas pomeron 
induced subprocesses result in particle production at midrapidity dominantly. 
Following this conjecture, we conclude that the relative importance of diffrac- 
tive to pomeron induced subprocesses decreases with increasing centrality of 
the collision. In this respect it would be helpful firstly, to increase experimental 
coverage in the target and/or the projectile fragmentation regions. This study 
has been proposed by |33| , |3^| . Secondly, it is important that experiments apply 
centrality cuts, in nucleus+nucleus as well as in p+nucleus and in elementary 
particle collisions. 
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3.1 Nucleus nucleus reactions 

In the following 4 subsections we analyse data from nucleus nucleus collisions 
at 1. RHIC, 2. SPS, 3. AGS and 4. SIS energies. 



3.1.1 Central Au+Au collisions at ^=130 GeV 



We used 3 measured ratios at midrapidity from |14]]T5H : A/A, p/p and h~ / 
Ncharged and imposed strangeness conservation. We exclude weak decay prod- 
ucts as is done in the experiment [I4"|. The data from reference JTH] are pre- 



liminary. The predicted and the experimental ratios are shown in table [l|. The 
resulting x 2 /DOF is (1.41/1) (CL 23%) for a temperature of 168 ± 40 MeV. 
When adding a systematic error of 15% linearly to account for the experimental 
feeding uncertainties and the fact that the ratios are measured at midrapidity 
only, the resulting \ 2 /DOF is 0.27/1 (CL 60%). The errors on the temper- 
ature and X s are not well estimated in this analysis due to the insensitivity 
of the ratios used, to the temperature at fixed fugacities. We estimate the 
error on the temperature searching for a variation of the \ 2 /DOF by 1, while 
changing the fugacities and fixing each time the T by imposing strangeness 
conservation. In order to improve the fit quality the addition of other ratios 
more sensitive to temperature variations is needed, like K/tt etc., which will 
soon be measured by the experiments at RHIC. 

After defining the (T, fi s ) values describing the particle ratios produced in 
central Au+Au collisions at -^5=130 GeV at the chemical freeze-out we ex- 
trapolate to the T at zero fugacities (table ^). Because of the very small value 
of the fugacities, the extrapolated values do not differ much from the ones at 
finite fugacitites. 



The x /DOF of all reactions are shown in table |21 



ratio 


model 


data 


A/A 


0.718 


0.70 ± 0.25 


p/p 


0.662 


0.650 ± 0.090 


h-/Nch 


0.491 


0.432 ± 0.0504 



Table 1: RHIC Au+Au at y/s= 130 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 
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Pb 

GeV 


Ps T 
GcV GeV 


K Pe 

GeV/fm 3 


T 

GeV 


\ s (eq,Pe) 


0.0390 


0.0084 0.168 


0.443 0.401 


0.169 


0.422 




± 0.040 




± 0.040 


+0.105 -0.156 






Pn 

1/fm 3 


T 

GeV 


\ s (eq,p n ) 






0.453 


0.168 
± 0.040 


0.443 

+0.11 - 0.164 






Ps 

1/fm 3 


T 

GeV 


X s (eq,ps) 






2.83 


0.168 
± 0.040 


0.443 

+0.11 - 0.164 


Table 2: 


RHIC Au+Au at 


Js= 130 GeV. 







Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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3.1.2 Central nucleus nucleus collisions at ~ 200 A GeV 

Here we show results from central Pb+Pb collisions at 158 A GeV, and from 
S+S and S+A collisions at 200 A GeV. 

Central Pb+Pb collisions at 158 A GeV 



For central Pb+Pb collisions at 158 A GeV, we use the particle ratios from 
table I in [IB| to compare with our model predictions. We don't consider a 



chemical potential for isospin | IB| , therefore we don't consider e.g. the 7r /ty~ 



ratio. We did not use the ratios in |16[ for which feeding correction is included. 
We did not use the d/d ratio, since there is strong experimental evidence that 
d and d are formed in the thermal hadronic freeze-out through coalescence 



T0| , |351| . Furthermore, the p may not freeze out in the chemical hadronic freeze- 
out due to the high cross section for p + p annihilation to e.g. tc as noted in 



We introduce a systematic error of 14% (quadratically added) in the ratios 
including WA97 data, because they are measured at midrapidity only. We 
add also quadratically) a 10% systematic error to ratios including h~ due to 
uncertainty in the feeding correction. The predicted ratios are shown in table 
[3] together with the experimental data. 

The x 2 & t which strangeness is conserved is 16.1 over 12 degrees of freedom 
(CL=20 % x 2 /DOF=1.3A) and corresponds to T=162 +10 -28 MeV. The 
asymmetry of the error (defined as the T deviation 1 unit of \ 2 /DOF away) is 
due to the fact that the 16.1 is not at the minimum. The minimum of the \ 2 is 
11.9 over 11 DOF (CL=50%, x 2 /DOF=1.08) at T=154 +14 -16 MeV. There 
is a probable physics process which may induce an imbalance of strangeness in 
experiment, if the measurements do not include the target and beam fragmen- 
tation regions: in particular particles with negative S especially hyperons may 
be produced near beam rapidity. In ref. |]37| it is conjectured that strangeness 



is not balanced (with missing particles with S negative). However we will 
demand in the following exact strangeness conservation. 

After defining the (T, fi s ) values describing the particle ratios produced in 
central Pb+Pb collisions at 158 A GeV at the chemical freeze-out we extra- 
polate to T at zero fugacities. The resulting equivalent temperatures for equal 
entropy density (p s ), energy density (p e ), and number density (p„) are shown 
in table g. 
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ratio 


\ \ \ s-\ f~\ 1 

IIlOQei 


data 


{P P)/ 11 


n 1 

U. 1UO 




n/ir° 
Tj/71 


7 9 




K°/n- 


0.142 


0.125 ± 1.90E-2 


K° s /h~ 


0.124 


0.123 ± 2.64E-2 


A/h~ 


1.04E-1 


7.70E-2 ± 1.54E-2 


A/tfJ 


0.835 


0.630 ± 1.02E-1 


K + /K~ 


1.78 


1.85 ± 9.00E-2 


K+/K~ 


1.78 


1.80 ± 1.00E-1 


-+/A 


0.170 


0.188 ± 3.90E-2 


(S- + S+)/(A + X) 


0.115 


0.130 ± 3.00E-2 




0.193 


0.232 ± 3.30E-2 




0.193 


0.247 ± 4.30E-2 




0.357 


0.383 ± 8.10E-2 


n/E 


0.153 


0.219 ± 5.00E-2 



Table 3: SPS Pb+Pb at y/s— 17 GeV, most central events. 

Predicted versus experimental particle ratios for the best fit of our model. 
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Pb 

GeV 


Ps T 
GeV GeV 




Pe 

GeV/fm 3 


T 

GeV 


\ s (eq,p e ) 


0.239 


0.0517 0.162 
+0.010 
-0.028 


0.565 


0.408 


0.169 

+0.008 

-0.030 


0.422 

+0.025 

-0.118 








Pn 

1/fm 3 


T 

GeV 


A s (eg,p n ) 








0.444 


0.167 

+0.008 

-0.030 


0.418 

+0.025 

-0.117 








Ps 

1/fm 3 


T 

GeV 


\ s (eq,p s ) 








2.8049 


0.168 

+0.009 

-0.031 


0.420 

+0.026 

-0.116 



Table 4: SPS Pb+Pb at y/s= 17 GeV, most central events. 
Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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Central S+A collisions at 200 A GeV 



In table | we give the same quantities for central S+A collisions at 200 A 
GeV with A=Au,W,Pb. The data used are taken from reference |T7|]. Here 
we mainly used the p/p and K + / K~ ratios to find (T, /it,, p s ) while we also 
suppose exact strangeness conservation. The resulted predicted ratios are p/p 
= 0.131 and K + / K~ = 1.59 to be compared with the experimental data p/p = 
0.12 ± 0.02 and K + /K~ = 1.59 ± 0.15. The resulting temperature is T=166 
+ 10 -28 MeV. The errors are taken to be percentually the same as in the 
Pb+Pb system. We do not perform a full fit, as the resulting temperature 



is within the errors in agreement with reference [1I71 where a full fit to many 



ratios is performed. The calculated values in |17j are T=165 ± 5 MeV, /Xb=175 
± 5 MeV, /i s =42.5 ± 4.5 MeV. 



Pb 


Ps T 


A, 


Pe 


T 


A s (eg,p e ) 


GeV 


GeV GeV 




GeV/fm 3 


GeV 




0.191 


0.0426 0.166 


0.544 


0.424 


0.169 


0.422 




+0.010 






+0.08 -0.030 


+0.025 -0.118 




-0.028 






+0.08 -0.030 


+0.025 -0.118 






A s 


Pn 


T 


X s (eq,p n ) 








1/fm 3 


GeV 










0.464 


0.169 


0.422 










+0.08 -0.030 


+0.025 -0.118 








Ps 


T 


A s (eg,p s ) 








1/fm 3 


GeV 










2.92 


0.169 


0.422 










+0.08 -0.030 


+0.025 -0.118 


Table 5: 


SPS S+A at v 


/ s= 19 GeV, most central events. 





Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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Central S+S collisions at 200 A GeV 



In table |6| the results for the central S+S collisions at 200 A GeV are shown, 
data are taken from references [18|]. We use K®/A, (B — B)/h~, A/A, K®/h~ 
and K + / K~ as well as strangeness conservation. The resulting x 2 /DOF 
=1.95/3=0.65 (CL ~ 57%) imposing strangeness conservation is at T=182 
+ 19 -8 MeV. The minimum of x 2 (x 2 /^ , OF=1.6/3=0.53 CL 66%) is at T=186 
± 12 MeV. 



P'b 




T X s 


Op 


T 


A,feo, pp) 


GeV 


GeV 


GeV 


GeV/fm 3 


GeV 




0.218 


0.058 


0.182 0.654 


0.182 


0.187 


0.473 






+0.019 +0.065 




+0.021 


+0.047 






-0.008 -0.003 




-0.009 


-0.023 








Pn 


T 


X s (eq,p n ) 








1/fm 3 


GeV 










0.803 


0.186 


0.470 










+0.009 


+0.047 










-0.020 


-0.023 








Ps 


T 


\ s (eq,p s ) 








1/fm 3 


GeV 










4.98 


0.187 


0.473 










+0.021 


+0.048 










-0.009 


-0.023 



Table 6: SPS S+S at \^s= 19 GeV, most central events. 
Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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3.1.3 Central nucleus nucleus collisions at ~ 10 A GeV 

Here we show results from central Si+Au collisions at 14.6 A GeV, and from 
Au+Au collisions at 11.6 A GeV. 

Central Si+Au collisions at 14.6 A GeV 

In table [7| the results for central Si+Au collisions at 14.6 A GeV are shown. 
We use the ratios K/n, K + /K~, A/ A and <f)/ir from [Q as well as strangeness 
conservation. We did not use p in view of the large annihilation expected at low 
energy. The resulting x 2 /DO F =6. 8/ '2=3. 4 (CL ~ 3.3%) imposing strangeness 
conservation is at T=120 +8 -35 MeV. We note however that the best x 2 
ignoring strangeness conservation 

(X 2 / DO F=0.092/2 with CL 63%) is at ^ 6 /T=4.05, fi e /T=0.9 and T=102 
MeV. 



Pb 

GeV 


Ps 

GeV 


T 

GeV 




Pe 

GeV/fm 3 


T 

GeV 


K(eq,pe) 


0.540 


0.107 


0.120 

+0.008 

-0.035 


0.548 


0.170 


147.4 

+0.013 

-0.056 


0.347 

+0.052 

-0.250 


Pb 


Ps 


T 




Pn 

1/fm 3 


T 

GeV 


X s (eq,p n ) 










0.182 


0.141 

+0.013 

-0.053 


0.322 

+0.048 

-0.23 


Pb 




T 


K 


Ps 

1/fm 3 


T 

GeV 


Xs(eq,p s ) 










1.161 


0.143 

+0.013 

-0.054 


0.330 

+0.049 

-0.237 



Table 7: AGS Si+Au at \/s= 5.4 GeV, most central events. 
Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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Central Au+Au collisions at 11.6 A GeV 

The results for central Au+Au collisions at 11.6 A GeV are shown in tables 
||, |9[ We used K + / K~ , K/n, p/n + and K/A ratios from 0. The resulting 
X 2 /DOF=0. 99/2=0. 496 (CL ~ 60%) imposing strangeness conservation is at 
T=96 +4 -5 MeV. The minimum of x 2 (x 2 /DOF=0.92/2=0AQ CL 63%) is at 
T=95 MeV. 



ratio 


model 


data 


< K > /A 


0.562 


0.675 ± 0.144 


K+/K- 


5.65 


6.303 ± 2.55 


p/n + 


1.088 


1.098 ± 0.127 


< K > / < IT > 


7.299E-2 


8.080E-2 ± 1.420E-2 



Table 8: AGS Au+Au at \fs= 4.9 GeV, most central events. 

Predicted versus experimental particle ratios for the best fit of our model. 
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Pb 

GeV 


Ps T 
GeV GeV 


Xs 


Pe 

GeV/fm 3 


T 

GeV 


\ s (eq,p e ) 


0.563 


0.084 0.096 

± 0.005 


0.281 


0.484E-01 


0.121 
± 0.08 


0.233 
± 0.037 








Pn 

1/fm 3 


T 

GeV 


\ s (eq,p n ) 








0.602E-01 


0.113 
± 0.079 


0.196 
± 0.037 








Ps 

1/fm 3 


T 

GeV 


A s (eg,p s ) 








0.387 


0.116 
± 0.08 


0.210 
± 0.037 


Table 9: 


AGS Au+Au at 


V~s= 


4.9 GeV, most central events. 



Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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3.1.4 Central nucleus nucleus collisions at ~ 2 A GeV 



Central Ni+Ni collisions at 1.9 A GeV 

The results for central Ni+Ni collisions at 1.8 A GeV are shown in tables [TU], 
|j~T| . We used the measured K + /K~ and K/tt ratios from []2l| , |22| (figure 4.4) 
and used additionally the K/A ratio deduced as A ~ A — A = 2(K + — K~) 
and imposed strangeness conservation. We add a systematic error of 15% for 
the feeding experimental uncertainty linearly. The resulting \ 2 / ' DO F=2. 16/1 
(CL ~ 10%) imposing strangeness conservation is at T=0.044 +0.0023 -0.0002 
GeV. The minimum of \ 2 (x 2 /DO F '=0.025/1 with CL better than 90%) is at 
T=0.0448 GeV. 



ratio 


model 


data 


< K > /A 


0.373 


0.275 ± 6.72E-2 


K+/K- 


20.90 


20.86 ± 6.88 


< K > / < IT > 


3.16E-3 


3.20E-3 ± 9.10E-4 



Table 10: GSI Ni+Ni at \^s= 2.6 GeV, most central events. 

Predicted versus experimental particle ratios for the best fit of our model. 



20 



Pb 

GeV 


Ps 

GeV 


T 

GeV 


Xs 


Pe 

GeV/fm 3 


T 

GeV 


\ s (eq,p e ) 


0.678 


0.102 


0.044 


8.93E-3 


0.987E-03 


0.602E-01 


1.55E-02 






± 0.0023 


± 0.00335 




± 0.0114 


+0.0435 
-5.55E-3 










Pn 

1/fm 3 


T 

GeV 


A s (eg,p„) 










0.161E-02 


0.504E-01 
± 0.0096 


5.16E-03 
+0.0144 
-1.84E-3 










Ps 

1/fm 3 


T 

± eq,ps 

GeV 


\ s (eq,p s ) 










1.18E-02 


5.32E-2 
± 0.0104 


7.38E-03 
+0.0207 
-2.64E-3 



Table 11: GSI Ni+Ni at y/s= 2.6 GeV, most central events. 
Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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3.2 Electron positron reactions 

The primary qq production in e + e~ collisions is not thermal due to the hard 
nature of the primary 7, Z couplings. However as long as the final particle 
multiplicity is much higher than 2, it is conceivable that this fact does not 
matter and the subsequent fragmentation of quark-antiquark jets into hadrons 
is thermal. We used the K/tt, p/n, n/p and A/n ratios. 
Our results for e + e~ collisions at ^=91 GeV are shown in table |12[ We used 
the data from reference U in particular the initially produced particles from 
the Pei model (table 1 in ||). We used as systematic error of the feeding cor- 
rection performed in this model, the difference between the Pei and the Jetset 
model (also in table 1 in H). The resulting minimum x 2 / DOF=2. 52/3=0. 84 
(CL ~ 47%) is at T=145 +39 -45 MeV and X s is 0.338 +0.127 -0.2. 
In reference |38| a temperature of 160.6 ±1.7 ± 3.1 is given for e + e~ collisions 
at 91 GeV, however with a not acceptable y 2 /DOF of 60.8/21, having a CL 
of 10~ 5 . In this reference strangeness has been weigted by a factor 7,, which is 
a free parameter of the fit. In reference |§ a temperature of T= 142.4 ± 0.018 
MeV and A«=0.295 ± 0.006 are extracted in e + e~ collisions at 91 GeV. Within 



our large errors we agree with both reference [38] and M. Our best value of 



the temperature is however nearer to the results of reference 



Pb 


Ps 


T 


Xs 


Pe 


Pn 


Ps 


GeV 


GeV 


GeV 




GeV/fm 3 


1/fm 3 


1/fm 3 


0. 


0. 


0.145 


0.338 


0.152 


0.208 


1.25 






+0.039 -0.045 


+0.127 -0.20 









Table 12: e + + e~ collisions at \fs= 91 GeV. 

Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. 
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3.3 Hadronic reactions 



3.3.1 Proton proton reactions 

Here we show results from p+p collisions at y/s=YI and 27 GeV. 
Proton proton reactions at ^fs=27 GeV 

For p+p collisions at ^fs=Tl GeV we used 9 measured ratios from reference 
H namely K/ir, p/p, n/p, rj/ir , A/AT S °, K+/K~, A/A, 0/vr, A ++ /p and 
imposed strangeness conservation. We add a systematic error of 15% linearly 
to account for the experimental feeding uncertainties. The predicted and the 
experimental ratios are shown in table O. We find a temperature of 128 +5 
-9 MeV, with a \ 2 1 X>OF=71/5=14.2, which has a CL of the order 1(T 14 , 
therefore the fit is not acceptable (table [L4]). 



In reference [23] the authors discuss p+p collisions at 4 different y/s. For 
a/s=27 GeV considered here, the obtained ^/DOF is 136.4/27 with a CL 
much less than 10 -5 . The obtained temperature in J23| is 169 ± 2.1 ± 3.4. 

The bad quality of the fits of p+p collisions at -\/i=27 GeV renders the re- 
sulting temperatures questionable. Nevertheless we point out that the two 
analyses yield incompatible temperature values within the errors. The analy- 
sis of P3| does not introduce systematic errors as we do. However, it allows 
for an arbitrary weighting factor (called 7 S ) acting only on strange particles. 
This factor is similar to an increase of the systematic experimental errors of 
strange particle multiplicities only. 
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ratio 


model 


data 


K/tt 


0.120 


7.35E-2 ± 1.80E-2 


p/p 


2.65E-2 


5.25E-2 ± 1.36E-2 


ir/p 


5.51 


3.14 ± 0.846 


eta/ n° 


6.89E-2 


7.75E-2 ± 3.15E-2 


A/K° 


0.533 


0.539 ± 0.158 


K + /K~ 


1.52 


1.48 ± 0.364 


A/A 


3.88E-2 


0.160 ± 6.75E-2 


0/vr 


1.07E-2 


5.034E-3 ± 1.34E-3 


A++/p 


0.148 


0.182 ± 4.97E-2 



Table 13: p+p collisions at y/s— 27 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. The 
X 2 /DOF of this fit is not acceptable. 
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Hb (GeV) 
GeV 


Us (GeV) 
GeV 


T (GeV) 
GeV 


A s 


Pe 

GeV/fm 3 


T 

GeV 


\ s (eq,p e ) 


0.243 


0.0282 


0.128 

+0.005 

-0.009 


0.329 


0.830E-1 


0.132 

+0.005 

-0.010 


0.283 

+0.0218 

-0.045 










Pn 

1/fm 3 


T 

GeV 


\ s (eq,p n ) 










0.123 


0.131 

+0.005 

-0.010 


0.279 

+0.0215 

-0.044 










Ps 

1/fm 3 


T 

GeV 


X s (eq,p s ) 










0.744 


0.132 

+0.005 

-0.010 


0.283 

+0.0218 

-0.045 



Table 14: p+p collisions at y/s= 27 GeV. 

Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. The x 2 'I DOF of this fit is not acceptable. 
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Proton proton reactions at y^i=17 GeV 



We used 6 measured ratios from references [0)111 namely K/tt, p/p, ff/p, 
K + /K~, <f>/ir, K/p and imposed strangeness conservation. We add a system- 
atic error of 15% linearly to account for the experimental feeding uncertainties. 
The resulting x 2 /DOF is (322/4)=80.6 (CL of the order 10" 66 !) for a tem- 
perature of 144 MeV, is completely unacceptable. 

The predicted ratios in table [15] are overestimating K + , K~ and <fi yields, 
therefore suggesting that if thermal conditions prevail, then at least two ther- 
mal reservoires are present. This may be due to the influence of diffractive 
processes in minimum bias triggers. 

We therefore repeat the fit with only ratios which are not much influenced by 
diffractive processes namely K/tc, K + /K~, and <p/ir. Note that the temper- 
ature at the minimum of y 2 does not satisfy strangeness conservation. This 
is combatible with the hypothesis of the importance of diffractive processes 
showing up in leading baryons (e.g. p, A). Imposing strangeness conservation 
the fit results in a temperature of 144 MeV with a bad x 2 /DOF=lQ9/l (CL 
~ 10~ 39 ). At the minimum of the \ 2 the temperatures is 96 +8 -9 MeV with 
a better x 2 /DOF=5.13/l (CL ~ 2%). These results are shown in tables [16 



T7j The very small errors resulting for p+p collisions at ^Js—Yl and 27 GeV 
should not be trusted because of the bad CL of the fit. 



ratio 


model 


data 


K/tt 


0.136 


5.97E-2 ± 9.90E-3 


p/p 


5.97E-2 


4.68E-2 ± 8.06E-3 


ix/p 


4.261 


3.59 ± 0.558 


K+/K- 


1.55 


1.61 ± 0.277 


0/vr 


1.83 


4.25E-3 ± 1.17E-3 


K/p 


0.578 


0.214 ± 3.56E-2 



Table 15: p+p collisions at a/s= 17 GeV. 

Predicted versus experimental particle ratios for the worst fit of our model. 



26 



ratio 


model 


data 


K/ir 


0.0723 


5.97E-2 ± 9.90E-3 


p/p 


5.36E-2 


4.68E-2 ± 8.06E-3 


7v/p 


47.7 


3.59 ± 0.558 


K+/K~ 


1.60 


1.61 ± 0.277 


0/vr 


2.06E-3 


4.25E-3 ± 1.17E-3 


K/p 


3.44 


0.214 ± 3.56E-2 



Table 16: p+p collisions at \fs= 17 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 
Only three ratios are fitted ((f>/ir, K + /K~, K/ir). 



Pb Ps 

GeV GeV 


T 

GeV 


K 


Pe 

GeV/fm 3 


T 

GeV 




\ s (eq,p e ) 


0.222 0.0277 


0.096 

+.008 

-.009 


0.151 


0.0124 


0.096 

+0.008 

-0.010 




0.121 

+0.034 

-0.039 








Pn 

1/fm 3 


T 

GeV 




\ s (eq,p n ) 








0.0293 


0.096 
+0.008 


-0.010 


0.121 

+0.034 -0.039 








Ps 

1/fm 3 


T 

± eq,p s 

GeV 




A s (eg,p s ) 








0.157 


0.097 
+0.008 


-0.010 


0.125 

+0.035 -0.04 


Table 17: p+p 


collisions at ^fs- 


= 17 GeV. 









Thermodynamic parameters for the best fit and temperatures and A s extra- 
polated to zero fugacities. Only three ratios are fitted ((f) /n, K + /K~, K/n). 
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3.3.2 Proton antiproton reactions 

Here we show results from p + p collisions at y / i=900 GeV and central and 
peripheral p + p collisions at y/s=1.8 TeV . 

Central proton antiproton reactions at ^=900 GeV 

For p + p collisions at 1/3=900 GeV we used 5 measured ratios out of the 
multiplicities from reference ||23|| . in particular K® /charged and N/ charged, 
A/ K^, A/Er, E~ / charged. The predicted and the experimental ratios are 
shown in table [18|. The temperature from the fit is: 

T(p + p v/i=900 GeV) = 143 +21 -37 MeV 

with a x 2 /DOF of 4.48/4 having a CL of ~ 34%, while A s is 

A s (p + P ^=900 GeV) = 0.330 +0.077-0.166. 

We add linearly a systematic error of 15% to account for experimental feed- 
ing uncertainties. The experimental errors are very large (e.g. 57% for E~) 
therefore the resulting temperature has also a large error. 



ratio 


model 


data 






K® /charged 


5.73E-2 


3.85E-2 


± 


1.06E-2 


n/ charged 


2.17E-2 


2.81E-2 


± 


1.25E-2 


A/K» 


0.181 


0.277 ± 


0. 


114 


A/E- 


7.96 


10.86 ± 


8. 


344 


E~ /charged 


1.305E-3 


9.83E-4 


± 


7.12E-4 



Table 18: p+p collisions at \fs= 900 GeV. 

Predicted versus experimental particle ratios for the best fit of our model. 
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Peripheral proton antiproton reactions at y / s=1.8 GeV 



For the most peripheral p + p collisions at ^=1.8 TeV we used 2 measured 
ratios from reference [^6| namely K/tt and p/n (table [19]). We use the only 2 
available ratios from experiment, which with zero fugacities and temperature 
as the only parameter leaves one degree of freedom. The most peripheral 
collisions are defined as the ones with the lowest measured charge multiplicity. 
The temperature from the fit is 

T(peripheral p + py/s=1.8 Te)= 140 ± 8 MeV 

with a x 2 /DOF of 18.04/1 having a CL of ~ 2.0E-4, therefore the fit is very 
bad. We add a systematic error of 0.004 linearly to account for the experi- 
mental feeding uncertainties. This error has been estimated from deviations 
between the ratios found in different experimental runs with p + p collisions 
at ^/s—1.8 TeV |5J|. The predicted and the experimental ratios are shown in 



table [19[ Then the resulting temperature is 154 ± 9 MeV and the x 2 / DOF is 
6.36/1 having a CL of ~ 1.1%. The quality of the fit is low as in the minimum 
bias p+p data at lower energies. X s is 

A s (peripheral p + p </s= 1.8 Te)= 0.372 +0.031 -0.034. 



ratio 


model 


data 


K/tt 


0.144 


0.115 ± 1.20E-2 


p/tt 


6.87E-2 


7.40E-2 ± 9.00E-3 



Table 19: Most peripheral p+p collisions at y/s— 1.8 TeV. 

Predicted versus experimental particle ratios for the best fit of our model. 
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Central proton antiproton reactions at ^=1.8 TeV 



The results for the most central p + p collisions at ^=1.8 TeV are shown in 



table The most central collisions are defined as the ones with the highest 
measured charge multiplicity. For the most central p + p collisions at y^i=1.8 
TeV we used again the 2 measured ratios from reference ]2B| namely Kj ir and 
p/n. The predicted and the experimental ratios are shown in table |2^. The 
temperature from the fit is 

T(central p + p y/s=1.8 TeV)= 150 ± 9 MeV 
with a x 2 /DOF of 2.28E-3/1 having a CL of ~ 96%. \ s is 

A s (central p + p ^=1.8 TeV)= 0.357 +0.035 -0.033. 
We did not add a systematic error. 



ratio 


model 


data 


K/tt 


0.1424 


0.1420 ± 2.00e-2 


p/n 


6.24e-2 


6.20e-2 ± 9.00E-3 



Table 20: Most central p+p collisions at \fs= 1.8 TeV. 

Predicted versus experimental particle ratios for the best fit of our model. 



We think that the dramatic change in the quality of the fit between peripheral 
and central p+p collisions can not be explained by overestimated experimental 
errors in the central collisions respectively by underestimated experimental 
errors in the peripheral collisions alone. Barring the obvious possibility that 
peripheral p + p collisions are non thermal the situation supports our previous 
conjecture on the important role of diffractive processes in peripheral and 
minimum bias hadronic collisions. 
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collision 


X 2 /DOF (CL) 


X 2 /DOF (CL) no SC 


Au+Au ^=130 GeV 


1.41/1 (CL 23%) 




with syst. error 15% 


0.27/1 (CL 60%) 




Pb+Pb v/i=17 GeV 


16.1/12=1.34 (CL 20%) 


11.9/11=1.08 (CL 50%) 


S+A v/i=19 GeV 


no fit performed 




S+S v^=19 GeV 


1.95/3=0.65 (CL 57%) 


1.6/3=0.53 (CL 66%) 


Si+Au ^=5.4 GeV 


6.8/2=3.4 (CL 3.3%) 


0.092/2=0.046 
(CL 63%) 


Au+Au v/s=4.9 GeV 


0.99/2=0.50 (CL 60%) 


0.92/2=0.46 (CL 63%) 


Ni+Ni Vs=2.6 GeV 


2.16/1 (CL 10%) 


0.025/1 (CL 90%) 


p+p ^=17 GeV 


169/1 (CL 10- 39 ) 


5.13/1 (CL 2%) 


min. bias (3 ratios) 






p+p ^i=27 GeV 


71/5=14.2 (CL 10~ 14 ) 


68.6/5=13.7 (CL lO" 13 ) 


min. bias 






p + p ^=900 GeV 


4.48/4 (CL 34%) 




min. bias 






p + p y/s—1.8 TeV 


2.3E-3/1 (CL 96%) 




p + p ^/s—1.8 TeV 






peripheral 


18/1 (CL 2E-4) 




no syst. error 






periph. 


6.36/1 (CL 1.1%) 




with syst. error 






e+ e - v^=91 GeV 


2.52/3=0.84 (CL 47%) 




minimum bias 







Table 21: Chisquare per degree of freedom and confidence level (CL) for each 
collision system analysed. The reactions are taken with central trigger unless 
differently stated. Strangeness conservation is taken into account in the second 
row while it is not in the 3d row (denoted by SC=strangeness conservation). 
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3.4 The initial energy density estimation 



methods. The first method (a) is based on Bjorken's formula [40 



To estimate the initial energy density achieved in each collision we use several 

(13) 



{dE T /dy) 

ycm 



where (dE T /dy) ycm is the transverse energy (E T = \Jp\ + m 2 ) per unit rapid- 
ity at midrapidity, Rt is the transverse radius of the particle source after a 
formation time r taken to be 1 fm/c. 

We estimate the (gL&t / 'dy) ycm when not available from experiment, as follows. 
We first estimate the total Et as 



(dE T /dy) tot (ycm) = ^ [ (dN/dy) a (ycm)yJ (pr) 2 a + m* ] (14) 

a 

(Method (a)) 

where the sum runs over particles of type a with multiplicity per rapidity 
unit at midrapidity (ycm) (dN/dy) a , mean transverse momentum (pt) q an d 
mass m a . The mean transverse momentum is also taken at midrapidity, when 
available. We take three types of particles: it, K and nucleons, as well as 
antinucleons when available. To calculate the differencial (dN/dy) ycm , when 
not available, we first find the mean (dN/dy), dividing the total particle mul- 
tiplicity by the total rapidity interval and then we use an extrapolation factor 
A to extract the value at midrapidity. This factor is found from measured ra- 
pidity distributions of particles at or as near as possible to the \/s considered. 
This may vary from nuclear to hadronic reaction and with y/s. In the cases 
where the dN/dy of the particles are measured at midrapidty we use these 
values. 

The formation time is taken in the literature usually in all reactions 1 fm/c, 
and this is what we use. In the following we give some examples of our e j n 
estimation with method (a). 

1. e + + e~ For the e + e~ collisions at \/s=91 GeV we used the multiplicities 
from ||. For two jet events in e + e~ collisions the jet axis is the longitudinal 
axis defining rapidity and pt refers to it. We used the dN/dy at midrapidity 
of ii, K, p in quark jets produced in e + e~ collisions at \/s=91 GeV by the 
DELPHI collaboration [||. 



32 



We used the mean transverse momenta from hadron hadron collisions at 
high y/s. In particular we use the mean transverse momenta from p + p colli- 
sions at 1.8 TeV from at charged multiplicity N~ 45 (to represent minimum 
bias values) namely 

(p T )„= 0.34 GeV, 

(pt) k — 0.5 GeV and 

(p T ) p = 0.59 GeV. 

In hadron collisions the mean transverse momenta do not change much with 
yfs above 10 GeV (see e.g. figure 2.7 in p2| ). The transverse size of the initial 



hadronic system of quarks and antiquarks produced, is given by the uncertainty 
principle, as 1 over the average px- For the transverse radius of e + e~ collisions 
we take therefore the inverse mean p? of 0.34 GeV, giving Rt=0.6 fm. The 
resulting initial energy density is: 

e in (e + e- y/s=91 GeV)=1.84 GeV//m 3 . 

The systematic error (which in this case is dominated by the uncertainty in 
the Rt determination) is about 50%. 

2. p + p at a/s=1.8 TeV For p + p we use the same method, while we 
use as transverse radius the radius of the nucleon of 0.8 fm. We use again 
pions, kaons, nucleons and antinucleons, and the mean transverse momenta 
measured in f26fl . We derive the dEx / dy(ycm from the measured charged 



total multiplicity per unit rapidity at midrapidity, using the ratios and mean 



transverse momenta measured as a function of charged multiplicity in p6| . No 
extrapolation factor to midrapidity is needed here. We find 

e in (dNjdr] = 19) = 7.61 GeV//m 3 , 

e in (dNjdr] = 15) = 5.82 GeV//m 3 , 

e ^ (dNjdr) = 11.5) = 4.31 GeV/fm 3 , 

e ^ (dNjdr] = 5.4) = 2.03 GeV//m 3 , 

e ^ [dNjdr] = 3) = 0.772 GeV//m 3 . 
In reference the e i n for dN c /di] = 15 of 3 GeV/fm 3 is calculated using 
pions only, and therefore is lower than our estimate^]. 

3. p + p ^=900 GeV For p + p at ^=900 GeV we use the multiplicities 



from []38| and the mean transverse momenta for p+p at ^/s—1. 8 TeV from 
2 We find the same e i n as reference p6[, when we use only pions. 
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We took the mean multiplicity per unit rapidity without extrapolation factor 
A to midrapidity since we don't know the shape of the rapidity distributions. 
We find 

e in (p + P Vs=900 GeV)= 1.23 GeV//m 3 

The pseudorapidity distribution of charged particles produced in p + p colli- 
sions at y/s=1.8 TeV is relatively flat over 6 rapidity units as measured by 
CDF ||49|| . However this may differ at the lower y/s of 900 GeV. We use the ex- 
trapolation factor A=1.49 (=dN/dy(ycm)/dN/dy(mean)) from the measured 
dN/dy distribution in ||l). The resulting e has an uncertainty of about 50%. 



4. Au+Au ^=130 GeV For Au+Au at ^=130 GeV we used the mean 
Pt fromp+p collisions at y/s=1.8 TeV and the particle multiplicities at midra- 
pidity (negative hadrons, total charged multiplicity, nucleons estimated from 
p — p ~ N charged — ^negatives-, as well as antinucleons using the p/p ratio) from 
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With the above method (a) (equations [13] , [14]) we calculated the e j n for 



several reactions. We summarize the results here: 

e in (Au+Au -\/i=4.9 GeV) = 0.86 GeV/ fm 3 , with mT — m Q instead of 
for nucleons 

e j n (Si+Au y/s=5A GeV) = 0.53 GeV/ fm 3 , with m T — m instead of ttit 
for nucleons 

In the following we used tut for nucleons too 

e in (Au+Au v^=4.9 GeV) = 1.70 GeV//m 3 , 

e m (Si+Au v^=5.4 GeV) = 1.53 GeV//m 3 , 

e m (S+S v/i=19 GeV) = 1.21 GeV//m 3 , 

e m (Pb+Pb v^=17 GeV) = 2.4 GeV//m 3 , 

e in (Au+Au v^=130 GeV) = 6.34 GeV//m 3 , 

e in{p + P Vs=900 GeV) = 1.23 GeV//m 3 , 

e in{p + P \/s=1.8 TeV most central) = 7.61 GeV/fm 3 , 

£ in{p + P \/s=1.8 TeV most peripheral) = 0.77 GeV/fm 3 , 

e m(p+P \/s=17 GeV) = 0.42 GeV//m 3 , 
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e m (e+ + e" v/s=91 GeV) = 1.84 GeV//m 3 



The initial energy density has been estimated by experimenters for some re- 
actions shown here, using an other method (/?), namely the Bjorken formula 



(equation 13) and the measured transverse energy 



E T (lab) = (Esin6) 



lab 



(15) 



(method 

where E is the total energy measured with e.g. calorimeters and 9 is the angle 
to the incident beam direction. The resulting values, (used in figures 1, 2 and 
5 in the following section), are: 

e m (Pb+Pb v/i=17 GeV) = 3.2 GeV/fm 3 pfQ, 

e in (S+Au) v/i=19 GeV) = 2.6 GeV//m 3 | 

e in (S+S v/i=19 GeV) = 1.3 GeV//m 3 ||,|4 



e in (Au+Au v^=4.9 GeV) = 1.3 GeV//m 3 g|, 

e j n (Si+Au -^=5.4 GeV) = 0.9 GeV/ fm 3 (see |0| and references therein). 



We estimate the maximal initial energy density with a third method (7) , taking 
the nuclear energy density of two overlapping nuclei 2 times the 7 factor of 
the colliding particles in the center of mass minus one: 



e 7 = 2 e A (7 - 1) (16) 

(method (7)) 

with 7 = (y/s 1 2) I m nuc i eon , and €A,smaii =0.1 79 GeV/fm 3 for small nuclei and 
e A,fejg=0.138 GeV/fm 3 for large nuclei is the normal nuclear matter density. 
The value in equation [16] multiplied by the stopping power gives an estimate 
of the initial energy density. 

We apply this method to the Si+Au and Au+Au collisions at \fs ~ 5 GeV, 
Ni+Ni collisions at \/s=2. 6 GeV, and Pb+Pb collisions at \/s=17 GeV, to 
calculate the maximal achieved initial energy density, yielding: 

3 The first three values here are valid for head-on collisions Esl. 
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e m (Au+Au v/i=4.9 GeV) = 0.44 GeV//m 3 , 
e m (Si+Au v^=5.4 GeV) = 0.67 GeV//m 3 , 
e m (Ni+Ni ^=2.6 GeV) = 0.276 GeV//m 3 , 
e in (Pb+Pb y/s=17 GeV) = 2.25 GeV//m 3 . 



Few comments: 



1. A problem with method (a) is, that in general it underestimates slightly 
the E t , since not all particles are taken into account. 

2. A general problem is with equation [15] of method ((3): For particles 
produced near midrapidity E x sinO is approximately pr not Et, for 
7cm >> 1, whereas E x sinO is approximately my for particles moving 
nonrelativistically in the lab frame. 

3. A problem with method (7) is, that it assumes geometrical and not 
dynamical compression, possibly slightly underestimating the e j„. 

4. The methods (a) and (3) should in principle give similar results. Our 
e j n (method (a)) agrees well with the calorimetric estimation of NA35 
for S+S at v^=19 GeV and of NA49 for Pb+Pb at y/s=17 GeV \MM 



for the 5% a to t centrality trigger. We compare to e in (Pb+Pb) ~ 0.77 x 
3.2 GeV//m 3 = 2.46 GeV//m 3 , with 0.77 = (dE T /drj)(h%a tot )/(dE T /drj) 
(head on collisions) from |43|. This correction is needed because the 
e i n =3.2 GeV//m 3 is estimated for head on collisions, while our calcula- 
tion uses m T and d/dy's taken with the 5% a tot centrality trigger. 
Additionally, also method (7), - and therefore all three methods (a), (/3) 
and (7)-, give similar results for Pb+Pb collisions at a/s=17 GeV. 



5. We find however different e values for the low -Js reactions Si+Au 
and Au+Au at 5.4 and 4.9 GeV, for each one of the above methods with 
a scatter of a factor of two or more. Also the maximal estimated e i n 
with method (7) is lower than the results from methods (a) and (j3). 
We therefore conclude that the Bjorken estimate may not be adequate 
for low energy reactions. This estimate was not meant to be used in the 
nonrelativistic regime. The method (7) seems more adequate for these 
systems. Nevertheless we show two distinct cases in the following figures. 

Using method (7) the e j„ is slightly higher for Si+Au at 5.4 GeV than for 
Au+Au at 4.9 GeV y/s, unlike the results of methods (a) and {(3). The 
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stopping power which is missing in the equation [16], can however hardly 
be very different for these two reactions, but may have an influence on 
the ordering of the e $ n in the above two reactions. 



Systematic error estimation summary Comparing results from methods 
(a) and (7) we estimate the systematic error in Si+Au and Au+Au collisions 
at ^=5.4 and 4.9 GeV to be about 50%, mainly resulting from the difficulty 
in applying the Bjorken formula. For the higher energy nucleus+nucleus colli- 
sions the systematic error is smaller ~ 30% [|ll|] resulting from the estimation 
of Et- The systematic error for e + + e~ is about 50%, resulting from un- 
certainties in the transverse area definition. The systematic error for p + p 
colllisions at y/s=900 GeV results from the Et definition and the extrapola- 
tion to midrapidity and is estimated to be about 50% . For p + p collisions at 
the Tevatron, the uncertainty comes mainly from not taking all particles into 
account with method (a) and is about 30%. The systematic error on the e j n 
estimate for Ni+Ni at y/s=2.6 GeV is 50%, resulting from comparison of our 



e iri value to model calculations 146 



3.5 The combined results: T, A s at zero fugacities 

Figures || and |3| show the resulting temperature and the X s factor extrapolated 
to zero fugacities along isentropic paths, as a function of the initial energy den- 
sity. The latter is taken from the experimentally measured transverse energy, 
when available, that is, for Si+Au ^=5.4 GeV, Au+Au a/s=4.9 GeV, S+S 
v/i=19 GeV, S+A y/s=19 GeV and Pb+Pb yfs=\7 GeV. For the remaining 
colliding systems, where no initial energy density estimation is available, we 
use our method (a) based on equation |TJ], exept for the Ni+Ni system where 
we use the method (7). The fact that the initial energy density estimation in 
figures ^ and was not done with the same method, introduces an additional 
systematic error on the e in scale. To reduce this uncertainty and show the 
systematic error of the initial energy density estimation we use in the following 



the e in estimated for 1.) high y/s (>10 GeV) by equations |TJ| and |TJ], and for 
2.) low yfi (<10 GeV) by equation 

Figures ^ and |5| show therefore the temperature and the X s factor again, ex- 
trapolated to zero fugacities along isentropic paths, as a function of the initial 
energy density, with a different estimation of the latter. In particular, the ini- 
tial energy density is not taken from the experimentally measured transverse 



energy, but is estimated using our method (a) based on equation [TJ], with the 
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exeption of the data at low y/s (< 10 GeV), that is, for Ni+Ni at a/s=2.6 GeV, 
Si+Au at v /s=5.4 GeV and Au+Au at t/s=4.9 GeV, where we used equation 



T6l which as discussed in the previous section, seems more adequate for the 



low energies as the Bjorken estimate. 

Figure [8] (a) shows the number density of kaons as a function of the initial 
energy density from reference |[L2|| , while in figure ^ (b) the temperature ex- 
trapolated to zero fugacities along isentropic paths, is shown as a function of 
the initial energy density. In figure || (a), the collision systems p+p at y/s=17 
GeV, S+S at js=l9 GeV, Pb+Pb at y/s=17 GeV, Au+Au at ^=4.9 GeV 
and Si+Au at -\/s=5.4 GeV are shown. In figure |8| (b) all the above reactions 
are shown and additionally the remaining analysed heavy ion collisions in the 
last section, namely S+A at y/s=19 GeV, Au+Au at -^5=130 GeV and Ni+Ni 
at y / i=2.6 GeV. The initial energy density of the common reactions displayed 
in figures [8] (a) and (b), -that is, of all reactions shown in |] (a)-, is defined in 
the same way, to allow their direct comparison. In particular, the initial energy 
densities for all heavy ion systems is taken from the experimental calorimetric 
measurements (method (3 above) when these are available. For Au+Au at 
•y/s=130 GeV we use our estimate with method (a). For p+p at y/s=17 GeV 



we used the estimate from reference [121 of 



e in (p+p at v/i=17 GeV) from @=0.85 GeV/fm 



This is found in [|l!| by a) estimating the dependence of the e j n for Pb+Pb 
collisions at ^/s—Yl GeV on the number of participant nucleons N and b) 
extrapolating this function to N=2. 

The rise and subsequent saturation seen in the kaon number density below 
and above e=1.3 GeV/fm 3 (figure |8| (a)), shows a clear relation to the same 
behaviour seen in the temperature (figure || (b)) as a function of the initial 
energy density. 

It is apparent from figures @, |3|, £|, ||, || and the present discussion, that a 
precise calculation of the critical energy density affecting all figures (0-§), needs 
more experimental data in the e region around 1 GeV/fm 3 . The errors on the 
temperature (figures 0, f|) are large around this region. An other uncertainty 
arises from the determination of the initial energy density especially at low y^i 
where the Bjorken estimate may not be adequate. Our best estimate of the 
initial energy density, is shown in figures [|, ||, || 0. In figures |3], |, |5|J^ we 
show results demanding a confidence level greater than 1%. In figures ^| and |7| 
we show the results demanding a confidence level greater than 10%. 
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Figure 2: Temperature at chemical freeze-out extrapolated to zero fugacities 
along an isentropic path, as a function of the initial energy density for several 
nucleus+nucleus, hadron+hadron and lepton+lepton collisions. We demand 
for the fits confidence level > 1%. 
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Figure 3: The number of strange antiquarks over the mean number of light 
ones at chemical freeze-out and for zero fugacities as a function of the initial 
energy density for several nucleus+nucleus, hadron+hadron and lepton+lepton 
collisions. We demand for the fits confidence level > 1%. 
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Figure 4: Temperature at chemical freeze-out and for zero fugacities as a func- 
tion of the initial energy density for several nucleus+nucleus, hadron+hadron 
and lepton4-lepton collisions. The initial energy density has been estimated in 
a different way as the one in figures ^| and |] (see text). We demand for the fits 
confidence level > 1%. 
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Figure 5: The number of strange antiquarks over the mean number of light 
ones at chemical freeze-out and for zero fugacities as a function of the initial 
energy density for several nucleus+nucleus, hadron+hadron and lepton+lepton 
collisions. The initial energy density has been estimated in a different way as 
the one in figures |2] and |3| (see text). We demand for the fits confidence level 
> 1%. 
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Figure 6: Temperature at chemical freeze-out and for zero fugacities as a func- 
tion of the initial energy density for several nucleus+nucleus, hadron+hadron 
and lepton+lepton collisions. The initial energy density has been estimated in 
a different way as the one in figures ^| and |] (see text). We demand for the fits 
confidence level > 10%. 
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Figure 7: The number of strange antiquarks over the mean number of light 
ones at chemical freeze-out and for zero fugacities as a function of the initial 
energy density for several nucleus+nucleus, hadron+hadron and lepton+lepton 
collisions. The initial energy density has been estimated in a different way as 
the one in figures |2] and |3| (see text). We demand for the fits confidence level 
> 10%. 
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Figure 8: (a) Kaon yield per interaction over the source volume at the thermal 



freeze-out as a function of the initial energy density (e j„) from reference [11,12 



(b) Temperature at chemical freeze-out and for zero fugacities as a function 
of the initial energy density for several nucleus+nucleus, hadron+hadron and 
lepton+lepton collisions. All collision systems shown in figure (a), are dis- 
played in figure (b) using the same initial energy density as defined in [n],[CJ]. 
We demand for the fits confidence level > 1%. 
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4 Conclusions 



We have discussed in section 2 the dependence of the critical temperature 
in the QCD phase transition on the vacuum pressure including in addition 
modes of noninteracting hadron resonances for u,d,s flavors of light quarks as 
defining the hadronic phase. Equating the pressure of this hadronic phase to 
the pressure of the plasma phase, represented by noninteracting u,d,s quarks 
and antiquarks and gluons as shown in figure 1, we obtain for zero fugacities 
T cr = 194 ± 18 MeV . 

We have performed a thermal analysis of yields in multiparticle production for 
13 reactions summarized in table |2~I| and discussed in detail in section 3. The 



intensive thermal parameters ( temperature and baryon as well as strangeness 
fugacity ) together with an error estimate are used to extrapolate the state 
associated with the chemical freeze-out of each reaction, studied along curves 
of equal entropy- , energy- and number-density, to zero fugacities. 
We represent the so obtained states by two intensive parameters 
T and A s = ..y,-^ , 

i.e. temperature and the ratio of antistrange to nonstrange (valence) antiquark 
abundances. These two quantities as functions of the initial hadronic energy 
density achieved in each reactions are displayed in figures @ - 0- The initial 
energy densities are estimated in subsection 3.4. The resulting error (horizontal 
error in figures |2] - is of the order of 50% around e crit , for the e + +e" collisions 
at y / i=91 GeV and for p + p collisions at ^=900 GeV, and approximately 
30% at higher e in . 

The reactions with an estimated confidence level above 10% were retained, 
which fall within the errors into two groups : 

group I with e 

in ^ 

lGeV/fm 3 

1. central Au+Au collisions at RHIC \/s =130/1 GeV 

2. central Pb+Pb collisions at \Js = 17 A GeV 

3. central S+A (A=Pb,Au,W) collisions ^ = 19 A GeV 

4. central S+S collisions at y/s = 19 A GeV 

5. e++ e" collisions at LEP ^ =91.19 GeV 

6. p + p collisions at \fs = 900 GeV 

7. central p + p collisions at \/s = 1.8 TeV 

group II with e 

in ~ \GeV I ffn 3 
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8. central Au+Au collisions at y/s = 4.9 GeV 

9. Ni+Ni collisions at GSI = 2.8 A GeV 



Both quantities T and A s saturate to constant values 



T lim = 155 ± 6 ± 20 MeV (syst) ; X Um = 0.365 ± 0.033 ± 0.07 (syst) 

(17) 

above the dividing energy density e j„ ~ 1 GeV/ fm 3 . Within the errors this 
is compatible with the critical energy density e cr u ~ \GeV j fm 3 obtained in 
lattice QCD j47j as well as with the results on the critical parameters obtained 



in section 2 ( eq. 12 ). This latter agreement and the extension of the work in 



|fi3| to include excited hadronic degrees of freedom is a new result. 



This limiting temperature is expected to be somewhat below the critical one, 
because hadronization is not an instantaneous process. The mean values show 
a difference in temperature of 40 MeV but the errors mainly due to the ap- 
proximations of free quasiexcitations and experimental errors are large adding 
linearly also amount to ~ 40 MeV. 

We estimate the systematic error on the limiting T, X s values, fitting the high 
e i n region of the figures f|, |5] using once a line fit and once a horizontal line fit, 
varying in both cases the number of bins taken, and estimating the deviations 
of the mean values resulting from each of the two fits. The error is ~ 20 MeV 
for the temperature and 0.07 for the X s . 

We interpret this saturation as characteristic for all reactions which before 
chemical freeze-out have passed through the quark gluon plasma phase, as 
opposed to those which remained throughout in the hadronic phase. 

We note that due to the large errors the dividing line as drawn at e „ ~ 
lGeV/ fm 3 is subject to a corresponding uncertainty. 

The saturation phenomenon observed here is in our interpretation confirmed 
by the same phenomenon derived for the kaon number density as a function 
of the initial energy density by one of us pT|JT2]| (figure ||). 
The agreement of the theoretical thermodynamic model with experimental 
ratios of identified particle multiplicities is quantified in the errors as discussed 
in section 3. We emphasize that agreement cannot be perfect due to several 
obvious and less obvious facts : 

a ) The rapidity and transverse momentum distributions do not agree with the 
spherically symmetric distributions in the thermodynamic description. 
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b ) There may well exist kinematically distinct thermal systems, as is indi- 



cated by the phenomenon of longitudinal and transverse flow | 31| 



The target and projectile diffractive regions may form distinct further 
thermal subsystems characterized by a different temperature than the 
particles produced near midrapidity (48 . 



The results derived here do indicate in our interpretation a high degree of 
thermalization of the systems studied. This implies the clear indication, that 
the quark gluon plasma phase is part of the hadronization process characteristic 
for hadronic, leptonic and nuclear reactions above the critical energy density. 
The universal excitation of the quark gluon plasma phase in nuclear and 
hadronic as well as leptonic multiparticle production is a new result. 

The behaviour of strangeness production as reflected by the parameter A s also 
reveals a new aspect: the saturation phenomenon corresponds to a strangeness 
enhancement relative to states with low initial energy density. No 'strangeness 
enhancement' is seen for heavy ion collisions relative to hadronic and leptonic 
reactions, when compared at the same initial energy density and zero fugacities. 
This enhancement arises only if systems with very different thermal properties 
in particular different fugacities are compared to one another at the same 
v^i, therefore at different energy density when nucleus nucleus collisions are 
compared to elementary particle collisions. 

LHC results for central p+p and Au+Au collisions at y/s— 14, 5.5 TeV will 

be important for the confirmation of the above picture. 

We conclude by proposing further experimentation in the following areas : 

i ) study of p+ p collisions at the Tevatron, with centrality selection on trans- 

verse energy, similar to nuclear collisions, and possibly with an improved 
particle identification option including the fragmentation regions. 

ii ) in conjunction with i) an extension to compare open and closed c+ c 

production. 

iii) within the program of heavy ion collisions at RHIC and at the SPS, 
allready under way, a dedicated study in the neighborhood of e injCrit = 1 
GeV/fm 3 , e.g. for e in =(0.6-1.5) GeV/fm 3 , in nucleus+nucleus, p+A and 
p+p collisions, using centrality selection. 
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